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We study production and equilibration of pions and muons in relativistic electron-positron-photon 
plasma at a temperature T <C m M , m^. We argue that the observation of pions and muons can be 
a diagnostic tool in the study of the initial properties of such a plasma formed by means of strong 
laser fields. Conversely, properties of muons and pions in thermal environment become accessible 
to precise experimental study. 
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I. INTRODUCTION 



The formation of a relativistic (temperature T in MeV range), electron-positron-photon e~,e + ,7 plasma (EP 3 ) 
in the laboratory using ultra-short pulse lasers is one of the topics of current interest and forthcoming experimental 
effort [l|, 0. The elementary properties of EP 3 have recently been reported, see Q, where typical properties are 
explicitly presented for T = 10 MeV. One of the challenges facing a study of EP 3 will be the understanding of the 
fundamental mechanisms leading to its formation. We propose here as a probe the production of heavy particles with 
mass m 3> T. Clearly, these processes occur during the history of the event at the highest available temperature, and 
thus information about the early stages of the plasma, and even pre-equilibrium state should become accessible in 
this way. 

We focus our attention on the strongly interacting pions tt^jTt (m^c 2 140 MeV), and muons /Lt ± (m /J c 2 ^ 106 
MeV), (in the following we use units in which k = c = h = 1 and thus we omit these symbols from all equations. 
Both, the particle mass, and plasma temperature, is thus given in the energy unit MeV.) These very heavy, compared 
to the electron (m e c 2 = 0.511 MeV), particles are as noted natural 'deep' diagnostic tools of the EP 3 drop. Of special 
interest is the neutral pion ir° which is, among all other heavy particles, most copiously produced for T <C m. The 
7T° yield and spectrum will be therefore of great interest in the study of the EP 3 properties. Conversely, the study of 
the in- medium pion mass splitting Am = m 7T ± — m^o — 4.594 MeV at a temperature T ^ Am will contribute to the 
better understanding of this relatively large mass splitting between 7r°and ir^, Am/m = 3.34%, believed to originate 
in the isospin symmetry breaking electromagnetic radiative corrections. 

However, given its very short natural lifespan: 

7r°^ 7 + 7, t° = (8.4±0.6)10~ 17 s. 

7T° is also the particle most difficult to experimentally study among those we consider: its decay products reach the 
detection system nearly at the same time as the electromagnetic energy pulse of the decaying plasma fireball, which 
is likely to 'blind' the detectors. 

This plasma drop we consider is a thousand times hotter than the center of the sun. This implies presence of 
the corresponding high particle density n, energy density e and pressure P. These quantities in the plasma can be 
evaluated using the relativistic expressions: 

Hi = J 9ifi(p)d 3 p, (1) 
e= [j29iEifi(p)d 3 p, E i = ^Jm 2 i +p 2 (2) 

i 

p = \j Eft(^-^V<(p)* 8 ' ( 3 ) 

where subscript 167, e~,e + , 7r°,7r + ,7r~, fi(p) is the momentum distribution of the particle i and and gi 

its degeneracy, for % = e~, e + , 7, /i + we have gi = 2, and gi = 1 for 7r°,7r~7r + . For a QED plasma which lives 
long enough so that electrons, positrons are in thermal and chemical equilibrium with photons, ignoring small QED 
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interaction effects, we use Fermi and Bose momentum distribution, respectively: 

•^ e± = e (u-p e ±u c )/T + p A = e u- Pl /T _ l> ( 4 ) 

The invariant form comprises the Lorentz-scalar u-p e , a scalar product of the particle 4- momentum pf with the local 
4- vector of velocity w M . In absence of matter flow and in the rest (in the laboratory) frame we have 



p? = (Ei,pl). (5) 



When the electron chemical potential v e is small, irT 3> v e , the number of particles and antiparticles is the same, 
n e - — n e + . Physically, it means that the number of e + e~ pairs produced is dominating residual matter electron yield. 
This is the case for all laboratory experimental environments of interest here, in which T > 2 MeV is achieved. We 
thus will set v e = in the following. 

It is convenient to parametrize the electron, positron and photon e~,e + ,7 plasma properties in terms of the 
properties of the Stcphan-Boltzmann law for massless particles (photons), presenting the physical properties in terms 
of the effective degeneracy g(T) comprising the count of all particles present at a given temperature T: 

y = e = g(T)oT*, W = g'{T)aT\ a = ^. (6) 

For temperatures T <C m e we only have in this case truly massless photons and g(T) ~ g'(T) ~ 2 7 . Once temperature 
approaches and increases beyond m e we find g ~ g'(T) ~ 2 7 + (7/8)(2 e - + 2 e +) = 5.5 degrees of freedom. In principle 
these particles acquire additional in medium mass which reduces the degree of freedom count, but this effect is 
compensated by collective 'plasmon' modes, thus we proceed with naive counting of nearly free EP 3 components. The 
factor 7/8 expresses the difference in the evaluation of Eq. ([3]) for the momentum distribution of Fermions and Bosons 
Eq. Q, with Bosons providing the reference point at low T, where only massless photons are present. In passing, we 
note that in the early Universe, there would be further present the neutrino degrees of freedom, not considered here 
for the laboratory experiments, considering their weak coupling to matter. 

In figureUJwe present both g(T) and g'(T), as a function of temperature T in form of the energy density e normalized 
by crT 4 , and, respectively, the pressure P, normalized by <jT 4 /3 . The g(T) jumps more rapidly compared to g'(T), 
between the limiting case of a black body photon gas at T < 0.5 MeV (g — 2) and the case g = 5.5 for 7, e~,e + , 
since the energy density also contains the rest mass energy content of all particles present. The rise of the ratio at 
T > 15 MeV indicates the contribution of the excitation of muons and pions in equilibrated plasma. We note that 
the plasma produced pions (and muons) are in general not in chemical equilibrium. The distribution functions which 
maximize entropy content at given particle number and energy content are ||: 

A= T; \ w±)e «WT_i' ^T.-VWM' (7) 

where T„or^±\ and T M are particles fugacities. For Tj — » the quantum distributions shown in Eq. ([7]) turn into the 
classical Boltzmann distributions, with abundance prefactor T<. 

In the case of interest here, when T < to, it suffices to consider the Boltzmann limit of the quantum distributions 
Eq.®, that is to drop the 'one' in the denominator. Using the the Boltzmann momentum distribution and taking 
the non-relativistic limit we have: 

^ ee n w = T.ijT^K/T) - T. (^) 3/2 e-^' T + ..., (8) 

where K 2 (and further below also K\) are the modified Bessel functions of integer order '2' (and '1' respectively). 

The particle densities are shown on right in figure [TJ The top solid line is the sum of n e + +n e - , which is marginally 
bigger than the photon density (dashed, blue) which follows below. We also include in the figure the sum density of 
muons +n^- (red, dashed), and the density of the neutral pion 7r° (bottom solid line). The chemical equilibrium 
corresponds to T w o( ff ±) = T M = 1 is used in figureUJon right, since this is the maximum density that can be reached 
in the buildup of these particles, for a given temperature. Both heavy particle densities appear comparatively small in 
the temperature range of interest. However, in magnitude they rival the normal atomic density (~ 10 2 /nm 3 )already 
at T = 4 MeV, and 5 MeV, respectively. This high particle density in the chemically equilibrated plasma explains the 
relatively large collision and reaction rates we obtain in this work. In turn, this opens the question how such dense, 
chemically equilibrated EP 3 state can be formed - we observe that colliding two ultra intense circularly polarized and 



FIG. 1: On left: the ratios g = e/erT 4 and g = 3P/aT 4 as a function of temperature T; on right: the equilibrium densities 
of electrons (blue, solid line), photons (green, dash-dot line), muons (red, dashed line), pions (blue dotted line) as functions of 
temperature T. 



focused laser beams on a heavy thin metal foil(s) is the current line of approach. Initial simulations were performed [4|]. 
Many strategies can be envisaged aiming to deposit the laser pulse energy in the smallest possible spatial and temporal 
volume and this interesting and challenging topic will without doubt keep us and others busy in years to come. 

As it turns out, even a small drop of EP 3 plasma with a size scale of lnm iSjgiven the high particle density, opaque. 
The mean free paths U of particles 'i' are relatively short, at sub nano-scale [3j: 




Where the reference energy values (31.1 and 27.5 MeV) correspond to the mean particle energy at T = 10 MeV. 
Photons are subject to Compton scattering, and electrons and positrons to charged particle scattering. In fact 
these values of U are likely to be upper limits, since Bremsstrahlung type processes are believed to further increase 
opaqueness of the plasma 0. In our considerations plasma particles of energy above 70 MeV are of interest, since 
these are responsible for the production of heavy particles. We see that the mean free path of such particles has also 
nm scale magnitude. 

We note that a EP 3 drop of radius 2nm at T = 10 MeV contains 13 kJ energy. This is the expected energy content 
of a light pulse at ELI (European Light Infrastructure, in development) with a pulse length of about At — 10~ 14 s. 
For comparison, the maximum energy available in particle accelerators for at least 20, if not more, years will be in 
head on Pb-Pb central collisions at LHC (Large Hadron Collider) at CERN, in its LHC-ion collider mode, where 
per nucleon energy of about 3 TeV is reached. Thus the total energy available is 200 fiJ, of which about 10%-20% 
becomes thermalized. Thus ELI will have already an overall energy advantage of 10 9 , while in the LHI-ion case the 
great advantage are a) the natural localization of the energy at the length scale of 10 _5 nm, given that the energy is 
contained in colliding nuclei, and b) the high repetition rate of collisions. 

As a purely academic exercise, we note that should one find a way to 'focus' the energy in ELI to nuclear dimensions, 
and scaling the energy density with T A up from what is expected to be seen at CERN-LHC-ion (T < lGeVY, we 
exceed T = 150 GeV, the presumed electro-weak phase boundary. Such consideration lead the authors of Refs. |l|, 
to suggest that the electro-weak transition may be achieved at some future time using ultra-short laser pulses. 

Returning to present day physics, we are assuming here that T near and in MeV range is achievable in foreseeable 





(9) 
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future, and that much higher values are obtainable in presence of pulses with At < 10 _18 s, cAi < 0.3nm. Hence we 
consider production processes for tt° , 7r ± , ^ for T < 50 MeV. We study here all two body reactions in EP 3 which lead 
to formation of the particles of interest, excluding solely ej — ► en , and the related e~e + — > 777°. The presence of a 
significant (1.2%) fraction of ir° — > e + e~7 decays implies that these related two body processes could be important 
in our considerations. However, these reactions involve the 7r off-mass shell couping to two photons, which needs to 
be better understood before we can consider these reactions in our context. 

We also do not consider here the inverse three body reactions e + e~7 — ► 7r°, since there is no exponential gain in 
using n > 2 particles to overcome an energy threshold, here m^o. The independent probability of finding n particles 
with energy m T o /n each is the same for any value of n: 

PxP 2 ....P n oc ( e - m -°/ nT y = e-^. (10) 

This resolves the argument that more particles could overcome more easily the reaction barrier, n-body reactions 
with n > 2 are in fact suppressed in EP 3 by the weakness of the electromagnetic (EM) interaction, since adding an 
EM- interacting particle to the reactions process requires an EM- vertex with a = 1/137. Thus microscopic reactions in 
EP 3 involving n > 2 are suppressed by a factor 100 for each additional EM particle involved in the reaction. This does 
not mean that a collective/coherent process of heavy particle production by many particles is similarly suppressed: 
for example fast time varying electromagnetic fields provide through E ■ B a collective source of 7r° . We defer further 
study of this production mechanism which requires multi MeV -1 range oscillation to be present in EP 3 . 

In the following section, we introduce the master equation governing the production of pions and muons in plasma 
and formulate the invariant rates in terms of know physical reactions. In section HTT1 we obtain the numerical results 
for particles production rates and reactions relaxation times which we present as figures. In section IIVI we discuss 
these results further and consider their implications. 



II. PARTICLES PRODUCTION RATES 



A. 7T° production 

7T° in the QED plasma is produced predominantly in the thermal two photon fusion 0: 

7 + 7^ vr°. (11) 

Much less probable is the production of 7To in the reaction: 

e-+e+^n°. (12) 

These formation processes are the inverse of the decay process of no . The smallness of the electro- formation of ttq is 
characterized by the small branching ratio in ttq decay B = r ee /r 77 = 6.2 ± 0.510 -8 . Other decay processes involve 
more than two particles. ir° can also be formed by charged pions in charge exchange reactions. However, in EP 3 in the 
domain of T of interest we find that at first the neutral pions will be produced. These in turn produce charged pions. 
Therefore we introduce the pion charge exchange process in the context of charged pion formation in the subsection 
III CI and since these can be important, we show these explicitely here as well. 

Omitting all sub-dominant processes, the resulting master equation for neutral pion number evolution is: 

1 dN„o d 4 W 77 _ >7r o d 4 M / 7r o_ 77 

(13) 



V dt dVdt dVdt 

T+7T~-*TT°+1T 



dVdt dVdt 

where N^o is total number of 7r°, V is volume of the system, (i 4 W 77 _ +7r n /dVdt is the (Lorentz) invariant ir° production 
rate per unit time and volume in photon fusion, and d^W^o^^^/dVdt is the invariant tt° decay rate per unit volume 
and time. Similarly, d 4 W 7r + n - _^ 7r o +7r o /dVdt is the pion charge exchange ir° production rate per unit time and volume 
while d 4 W 7T o +7T a^ 7T + 7T - /dVdt is the corresponding reverse reaction loss rate. 

We assume that in the laboratory frame the momentum distribution of produced ir° are characterized by the 
ambient temperature. Eq. ([SJ defines the relation of fugacity to the yield. This equation allows now to study the 
production dynamics as if we were dealing with a 7r° in a thermal bath, and to exploit the detailed balance between 
decay and production process in order to estimate the rate of tt° production. This theoretical consideration should 
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not be understood as assumption of equilibration of ir°, which could upon production escape from the small plasma 
drop. 

In [tJ the detailed balance relation is derived in detail, which takes the form 



_ x rf 4 W 7r o_ tyr _ T _ 2 d 4 W, 
dVdt ~ 7 dVdt 



This allows that Eq. (jT3|) can be written in the form 

1 dN„o 



V dt 



(T - T 7r o)_R 7r o — (T^o — T 7r± )i? 7r o w o < _ +7r + 7r - . (15) 



For T Wo — > T 2 , — > T^. ± = 1 we reach chemical equilibrium, the time variation of density due to production and decay 
vanishes. 

The charge exchange process rate (i? 7r o 7r a«_ >7r + 7r - , last in Eq. (JT5J)) balances the first contribution in Eq. (|44p . where 
it will be further discussed. The rate R n o can be written as 

d 3 p^ f d 3 p2 7 f d?pi 1 4 4 



K ° = J WT^rJ (2. Y2E 2 J (2. )32f; l7 (27r) * (Pi7+«.7"ftr)x 

J2 KP1.P2-, |Af|p.)| 2 /.(^)/ 7 (pi 7 )/ 7 (P2 7 )T; 2 T; 1 e tl ^/ T . (16) 

spin 

where for 7r° formation there was the factor (1 + f^) which we reduced using the relation 

lT/± = TrVWT /±; (17) 

where Fermi (/+) and Bose (/_) distributions are implied for particle i. Similarly, in the 7r°-decay case we replaced 
the two stimulated decay factors (1 + / 7 ) 2 in that way. Eq. (|16[) follows Including in Eq. (fT()|) the prefactors required 
by Eq. (TTJJ and recalling time reversal invariance, i.e. M = M'\ 

|(pi 7 P2 7 \M\p^)\ 2 = \fa \M\p ljP2j )\ 2 . (18) 

We realize that the result, Eq. (|16|) is manifestly symmetric for the two reaction directions. It is interesting to note 
that in Boltzmann limit all fugacities cancel in Eq. (fT6)) . 

We introduce the pion equilibration (relaxation) time constant by: 

dn^o/dr^a 

T^o = . (19) 

Note that when the volume does not change in time on scale of r^o (absence of expansion dilution) and thus T 
is constant, the left hand side of Eq.ffTS]) becomes dn^o/dt. Given the relaxation time definition Eq. (Ti"§|) the time 
evolution for of the pion fugacity for a system at fixed time independent temperature satisfies: 

T w°— 77- = X 7 - 1 tt« - (I ff o - I OT ±) . (20) 

at H^o 

When the charge exchange reaction can be ignored, for T 7r o(< = 0) = we find the analytical solution T OT a = 
T 2 (l — e~ t / r "°), justifying the proposed definition of the relaxation constant. 

We note that Eq. ptJ]) also describes the decay of a ir°. Therefore, up to small modifications introduced by the 
thermal medium (see discussion below), 

The 7T° production rate is thus related to the decay rate l/r° by the simple formula 

dn^/dT^a f m7I T\ 3/2 e' m -/ T 
R *° " r° " (21) 

where in the last expression we have used Eq. ([5]) in the limit m » T. It is important for the reader to remember 
that derivation of Eq (|2ip is based on detailed balance in thermally equilibrated plasma, and does not require chemical 
equilibrium to be established. 



T [MeV] 



FIG. 2: The ratios r^o /t° as functions of temperature T for relativistic Boltzmann limit (blue, dashed line) and for quantum 
distribution in chemical equilibrium, Y n = T 7 = 1 (green, solid line). 

Now we consider how and why r w o ~ r° . It turns out that there are both relativistic and quantum effects which 
contribute and they (nearly) cancel: the relativistic effect arises because T n o in Eq. (f2"Tj) is in lab frame while the known 
t^o is in the pion rest frame. In the relativistic Boltzmann limit the correction is obtained considering the related 
time dilation effect 0] is: 

T - T °° - T o K 2 (m^/T) 

^-TTh^-^Kiim^/Ty (22) 

where < 1/7 > is average inverse Lorentz factor. We find that this effect implies that r w o in the lab frame increases 
with temperature. This effect is shown by dashed (blue) line in figure [21 Furthermore, with increasing temperature 
quantum distribution functions for photons and for the produced particle need to be considered. This leads to the 
result shown as solid line (green) in figure [2] Thus in general r^o > r° , by up to 14%. 
We can further evaluate exactly the reaction rate Eq. (fT6|) [3| : 



(2^) 2 T° n 7o E n T; W r -l 



^ = 7^^l ^r^-C (23) 



where 



with 



*(*>=/ <V Y - le(o - fe0 - 1T ^ e ( a+ K)-i ' ™ 



a= y/K+A . fc * (25) 



2T 2T 



This integral for T 7 = 1 takes the form: 



,(b-a) _ -(a+b)\ 



*<P-) = ^T-T) I 6 + in I 1 + (1 _ e6 _ a) ] ] • (26) 



This exact result (blue, solid line) is compared to the approximate result Eq. (|2"Tj) (green, dashed line) in figure We 
note that it is hard to discern a difference on logarithmic scale, especially so for small temperatures where the only 
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FIG. 3: The 7i~o production rate (blue, solid line) and approximate rate from Eg, (|21[l (green dashed line) as functions of 
temperature T. 



(small) effect is the relativistic time dilation. This implies that it is appropriate to use the simple result Eq. (|21|) in 
the study of 7r° production. 

Before closing this section we note that we can use exactly the same method to extract from the partial width of 
the ttq — > e + e~ the reaction rate for the inverse process, which will be discussed below. All arguments carry through 
in identical and exact fashion replacing where appropriate the Bose by Fermi distributions and using Eqll71 



B. Muon production 

In the plasma under consideration, muons can be directly produced in the reactions: 

7 + 7 ^/i++//-, (27) 
e+ + e~ -> + (28) 

For reactions (|27[) and (I28D the master evolution equation developed for the study of thermal strangeness in heavy ion 
collisions applies [jj.l9lllCl.HH (compared to these references our definition is changed, their i?i2^34 — > i?i2^34/ (T1T2) 
in order to make the forward-backward symmetry explicit ) 

- (T„ - T^)i? 77WAI+A1 - + (Tg - T^)i? e+e -^ AI+M ~. (29) 



y fit 7 M f 1 • \~e ~ n> 

Like before for ir° we consider the master equation in order to find appropriate definition of the relaxation time 
constant for /j, production. In no way should this be understood to imply that muons are retained in the small 
plasma drop. In chemically equilibrated EP 3 the [i production relaxation time is defined by: 

r - - \m — dJh rt v <30) 

where a suitable choice is a = 1, 2 for t = 0, 00, respectively (see below). The form of Eq. (f30"|) assures that, omitting 
the volume expansion, i.e. the dilution effect, the evolution of the muon fugacity obeys the equation 

ar^ = l-T% T 7 = T e = l, (31) 
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which has for T M (t = 0) = the simple analytical solution 

Y M = tanh t/ar^. (32) 

For t — > oo, near to chemical equilibrium, T M — > 1 — e~ 2t / aT,i , while for i — > 0, at the onset of particle production 
with small T M we have = tj (ar M ) . Hence, near to chemical equilibrium it is appropriate to use a = 2 in definition 
of relaxation time Eq. (|30|) , while at the onset of particle production, more applicable to this work a more physical 
choice would be a = 1. However, following the convention, in the results presented below the value a = 2 is used. 
The invariant muon production rate in photon fusion as introduced above is: 

f d 3 Pll + f d 3 p M - f d 3 p ly f d 3 p 2l . .4 ,4/ \ 

+ "- = J wmrJ wm~J wn^-J wr^ {2n) 6 

^2 | (pijP2 7 |A^ T ^ M+M - I | 2 / 7 (pi 7 )/ 7 (p 27 )/ AI (p A1 +)/ Al (p Al -)T 7 " 2 T~ 2 e tl ' (p f ++:p A'- )/T (33) 

spin 

and the invariant muon production rate in electron-positron fusion is: 

p f d 3 p + /■ d 3 p - f d 3 p e + f d 3 p e - 4 4 , v 

^•--^ = y (2 , )32 ^ + y (27r) 3 2 ^_ y (27r)32i , e+ y 5 ^ -*v x 

XI |(Pe+Pe- |M e+e ^ M+M -|p M+ p M -)|Ve(p e ^/ e (Pe-)/M(PM+)/M(^-) T ^ 2 T; 2 e n ' ( ^ ++P -- )/T . (34) 
spin 

We note that in Eq. (|33[) and Eq. (j34f in the Boltzmann limit all fugacities cancel, and that the forward-backward 
reaction symmetry is explicit. Moreover, it is interesting to note that despite inclusion of quantum effects (Bose 
stimulated emission and/or Fermi blocking), when using rates as defined in this paper, we don't change the master 
population equation system arising for Boltzmann particles. The only modification is a slight fugacity dependence of 
rates presented in Eq. (16]), Eq. Q , Eq. (|34|) . 

The IMe+e-^+^j- | 2 differs from often considered heavy quark production \M q q^ c5 \ 2 [H, EH (m c » m q ) by 
color factor 2/9, and the coupling a s — > a of QCD has to be changed to a = 1/137 of QED. Then we obtain, based 
on above references: 

2^ l-Me+e-^M+M-l =5 e 87ra — , (35) 

where m = 106 MeV is the muon mass, electron and positron degeneracy g e — 2, and s, t, u are the usual Mandelstam 
variables: s = [pi + P2) 2 , t = (pz — Pi) 2 , u= (pj, — P2) 2 , s + t + u = 2m 2 . For the total averaged over initial states 
\M\ 2 for photon fusion we have 

I *Wm- I 2 = v> 2 * 2 (" 4 + 2 + 4 + -^-) + 4^ + 4^) - (36) 

' \ \m z — t m z — u J \m z — t m z — u ) m z — t m z — u I 

where degeneracy g 7 — 2. Near threshold s m 4m 2 , with t,u ~ — m 2 we find 

|M 77 ^ M+M - 1 2 = 64^ 2 « 2 , |M e+e -^ M+M - 1 2 = 32^ 2 a 2 . (37) 

The e + e~ — > reaction involves a single photon, and thus it is more constrained (by factor 2) compared to the 

photon fusion, which is governed by two Compton type Feynman diagrams. However, in the rate we compute below, 
the indistinguishability of the two photons introduces an additional factor 1/2, so that both reactions differ only by 
the difference in the quantum Bose and Fermi distributions. 

Integrals in Eq. (f3"3"|) and (f3~4"| can be evaluated in spherical coordinates. We define: 

1 1 

q = Pi+p 2 ; p = -(pi - p 2 ); q'=P3+Pi; p = 2<i?3 - Pi); (38) 

z-axis is chosen in the direction of ~~q — pi + p~2\ 



% = (9o,0,0,0), p^ = (p ,psm6,0,pcos9), p'^ = (p' ,p' sin 4> sin x,p' sin <j> cos x,p' cos <j>). 



9 



1 (47r)(27r) f°° f s ~ q ° 
i? e + e -( 77) ^+ M - = T— T (o Hifl / dq ° / d q I dp a I d P'o I d P I d P' I d(cos9) I d(cos< 



Now we obtain [l(J: 

~ 1 + / (2tt) 4 16 

f (,- (,§ + j)" 2 ) . (* - h - * + * (cos,- a („♦ - ^) 

^ |M e+e _ (77) ^ w | 2 T- 2 /^ (y +Po) //, (y -Po) T f T( 7 )/e( 7 ) (y +Po) /e( 7 ) (y ~ Po) e^p{q Q /T), (39) 



where q* = qy 1 j 1 . The integration over p, p', cos 0, cos can be done analytically considering the delta- functions. 

The other integrals can be evaluated numerically. For the case of indistinguishable colliding particles (two photons) 
there is additional factor 1/2 implemented by the value 1 = 1, while for distinguishable colliding particles (here 
electron and positron) 1 = 0. 



C. tt ± production 

tt^~ can be produced in ■kqttq charge exchange scattering: 

7T° +7T° -> 7T+ +7T-, (40) 

as well as in two photon, and in electron-positron fusion processes 

7 + 7^7T + +7T~, (41) 

e+ + e _ -> 7T+ + tt~. (42) 

We find that for production, the last two processes are much slower compared to the first, in case that 7r density 
is near chemical equilibrium. Similarly, the two photon fusion to two ir°: 

7 + 7^7r° + 7r , (43) 

turns out, as expected, to be much smaller than one n° production. It is a reaction of higher order in a and the 
energy is shared between two final particles. 

The time evolution equations for the number of are similar to Eq. ([29]): 

y Jj_ = (Yff - T^±)i? 7r o 7r o^ 7r + 7r - + (T 2 - T 2 r± )i? 77 ^ 7r + 7r - + (Tg - T 2 r± )i? e + e -^ 7r+7r - . (44) 
In order to evaluate the pion production rates in two body processes we use reaction cross section, and the relation flij l : 

Ri2~« + «- = S^T^T r dsa^^K^s/T), (45) 



32tt 4 1 + /.,., , \ ■ ' 

" a th v 

(compared to reference Q our definition is changed i?i2-+34 — > i?i2^34/(Ti Y 2 )) where 

M{s) = (s — (mi + m 2 ) 2 ){s - {mi - m 2 ) 2 ), (46) 

mi and m 2 , gi and g 2 , Ti and T 2 are masses, degeneracy and fugacities of initial interacting particles. 

For the respective three cross sections we use, all results valid in the common range s < 1 GeV 2 we consider here: 

• The cross section for charge exchange 7r°scattering reaction Eq. (|4T))) have been considered in depth recently [l5[ : 



16tt s-4M 



a = 



9 Vs- 4M 2 



^(4 0) -4 2) ) 2 ; (47) 



where dp ' — off = 0.27/M n ± This is the dominant process for charge pion production, subject to presence of 
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For process Eq. (|4ip . the cross section of production in photon fusion we use [l£ 

27ra 2 / 4m! x 1/2 



1 



(l/2s + mf0(V 4s + m v) 



(48) 



where my — 1400.0 MeV. As we will see from numerical calculations given the cross sections for 77 — > tt + tt 
resulting production rates will be smaller than the charge exchange 7r°7r° — > tt + it~ reaction. 



• For process Eq. 



the cross section of ^ production in electron - positron fusion we use [l7| : 

ira 2 (s - 4m 2 ) 3/2 



„5/2 



\F(s)\' 



The form factor F(s) can be written in the form: 
F(s) = 



mpTpd 



m 2 p — s + T p (m 2 /k^)[k 2 (h(s) 
where h'(s) = dh/ds and 



h(m 2 )) + k 2 h'(m 2 )(m 2 - s)] - im p (k/k p ) 3 r p {m p / 'a/s) : 



(49) 



(50) 



1 

4 S 



1/2 



1 2 



1/2 



2 k 
h(s) = -= In 

7T yjs 



2k 



2m T: 



trip = 775 MeV, T p — 130 MeV, d — 0.48. Given this cross section we also find that the rate of charged pion 
production is small when compared to 7r -charge exchange scattering. 



• For reaction 



we have [l8j: 
(7(77 -> 



7T 7T ) = 



o? \J s — Am 2 \ 
8ir 2 ^ ) 



JL f 

Js 



(7(77+ 7T- 7T U 7T U ), 



(51) 



where 



/ s = 2(ln 2 (z + /z_) - tt 2 ) + ^(ln 2 (z + /z_) + n 2 )' 



(52) 



and z± = (1/2)(1 ± v/s-4™ 2 ). 

The cross sections for 7r + 7r~ pair production, evaluated using Eqs. (|47|) . ([48]) and (j49l) are presented in figure [4] 
as functions of reaction energy y/s . Top solid line (blue) is for charged pions production in tt° scattering Eq. (|40| , 
the magnitude of this cross section being very large we reduce it in presentation by factor 1000; the dashed line is 
for 7T + 7r - production in photon fusion Eq. pdj) ; dash-doted line is for electron positron fusion Eq. (|42| . The bottom 
solid line (green) is for photon fusion into two neutral pions, Eq. (|51[) . The prediction for cr 77 _ +7r + w - is about 480 nb 
(data 420 nb) at the peak near threshold [Hj], which is in agreement with calculations presented here. The reaction 
(Eq. (|4"3"|) ) is much smaller than others and we do not consider this reaction further. We note that some of 
these results are currently under intense theoretical discussion as they relate to chiral symmetry. For our purposes 
the level of precision of here presented reaction cross sections is quite adequate. 



III. NUMERICAL RESULTS 



A. Particle production relaxation times 



In figure [5] we show relaxation time r for the different processes considered as function of temperature T £ [3,50] 
MeV. Because of the large difference in production rates which can be compensated by different densities of particles 
present (magnitudes of fugacities) we introduce partial relaxation time for each of the three reactions ir°ir° — > 7r + 7r~, 
77 — > 7r + 7r~ and e + + e~ — > tt + tt~: 

1 dn^± /dT T ± 1 dn 7T ± /dT n ± 

Z xt 7r 7r Q. { _> 7r + 7t — Z 



1 dn w ± /dT 7T ± 



(53) 
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FIG. 4: The cross section a for pion pair production, and pion charge exchange (solid top line), as functions of \fs < 1 GeV 2 . 

When T <C m, we can use the Boltzmann approximation to the particle distribution functions. Since in this limit the 
density is proportional to T the relaxation times doesn't depend on T. Moreover, even for T — ► 50 MeV, we have 
for muons e~ m / T ~ 1/3, thus quantum correlations in phase space remain small, and the Boltzmann limit can be 
employed. To account for small deviation from Boltzmann limit arising towards the upper limit of the temperature 
range we consider, that is at T ~ 50 MeV, we used the exact equations with Tj = 1 to calculate r for each case, value 
corresponding to the maximum density that can be reached for a given temperature, for which the quantum effect is 
largest. In addition to these three cases Eg. (fS"3")) we show in figure [5] the muon production relaxation time Eg. ([501 ). 
the two photon fusion into 7r° relaxation time Eg. (|19|) . a nearly horizontal line (turquoise, bottom), which is slightly 
greater than the free space 7r° decay rate. Finally, the thin dash-dot line at about 10 s times greater value of time is 
the electron-positron fusion into tt", Eq.(12). 

B. Rates of pion and muon formation 



In figure [5] we show on left as a solid (blue) line as a function of fireball temperature the rate per unit volume 
and time for the process 7 + 7 — * ir° , the dominant mechanism of pion production. The other solid line with dots 
corresponds to e + + e~ — * ir° reaction which in essence remains, in comparison, insignificant. Its importance follows 
from the fact that it provides the second most dominant path to ttq formation at lowest temperatures considered, and 
it operates even if and when photons are not confined to remain in the plasma drop. 

We improve the rate presentation on the right hand side in figure [SJ considering that the formation of a plasma 
state involves an experimentally given fireball energy content £ in Joules, we use Eq.® to eliminate the volume V 
at each temperature T: 



d 2 W 



1 d A W^ 



1 



dtd£ 



gaT 4 dVdt 



gaT 4 



For chemical nonequilibrium, replace g ~ * Tt,g(T). Considering the (good) approximate Eo. (|21[) we obtain: 



R' 1 



\2ttT) 



3/2 e -m„/T 



(54) 



(55) 
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FIG. 5: The relaxation time r for the different channels of pion and muon production (see box), as functions of plasma 
temperature T. 



We use units such that h = c = k = 1 and thus R' is a dimensionless expression. Recalling the value of these constants, 
the units we used for R' derive from MeV s=1.60310 -4 MJ fs. 

The other lines in figure [5] address the sum of formation rates of charged pion pairs (dashed, red) by all reactions 
considered in this work, ir° + n — > 7r + + 7r~,7 + 7 — > ir + + ir~ , e + + e~ — > 7r + + 7r~. We also present the sum 
of all reactions leading to either a charged pion pair, or muon pair (dot-dashed, green) lines, that is adding in . 
7 + 7 — > /i + + /i~,e + + e~ — > /i + + fi~ . The rationale for this presentation is that we do not care how a heavy 
particle is produced, as long as it can be observed. The dashed (red) line assumes that we specifically look for 
charged pions, and dot-dashed (green) line that we wait till charged pions decays, being interested in the total final 
muon yield. The tt° production rate (blue, solid line) is calculated using Eq.Q16p and yields on the logarithmic scale 
nearly indistinguishable result from the approximation Eq.([2"T]). For tt^ 1 production we refer to section Hi CI and for 
fj, production we refer to III Bl 

In table U we show the values of key reaction rates R and relaxation times r at T = 5 and 15 MeV. We note the 
extraordinarily fast rise of the rates with temperature, in some instances bridging 15 - 20 orders in magnitude when 
results for T = 5 and 15 MeV are compared. 

In order to understand the individual contributions to the different reactions entering the sum of rates presented 
above, we show as a function of temperature in the figure [7] the relative strength of muon pair (left) and charge pion 
(on right) electromagnetic (7 -I- 7, e + + e~) production, using as the reference the 7 + 7 — ► ir° reaction. The fi 
production rates are calculated using Eq.(|39p with \M\ 2 from Eq. (f3"5)) and Eq. fBl)]) respectively. This ratio is smaller 
than unity for T ^ 20 MeV. For larger T, the muon direct production rate becomes larger than tt° production rate. 
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FIG. 6: On left, the invariant pion production rates in units of nm 3 fs , as a function of temperature T. On right the 
production rate R' per Joule energy content in the fireball, in units of MJ _1 fs _1 , in both cases for reactions shown in the box. 



TABLE I: Values of rates, relaxation times for all reactions at T = 5 MeV and T = 15 MeV 





T = 5 MeV 


T = 5 MeV 


T = 15 MeV 


T 


= 15 MeV 


reaction 


r [as] 


i? [nm-'fe" 1 ] 


t [as] 


R 


[nm- 3 f S - 1 ] 


77 <-> 7T0 


88 


3.3 10 3 


95 




1.2 10 12 


e + e~ <-» fx fiT 


1.2 10 10 


3.2 10~ 3 


1.9 10 3 




1.5 10 11 


77 <-> n + n~ 


1.0 10 10 


3.7 10" 3 


1.3 10 3 




2.1 10 11 


7r°7r° «-» 7r + 7r~ 


2.9 10 12 


2.1 10" s 


4.6 10 2 




9.5 10 10 


77 <-> 7r + 7r _ 


6.4 10 13 


9.7 10" 10 


5.1 10 4 




8.710 s 




7.8 10 15 


7.9 10~ 12 


9.5 10 5 




4.6 10 7 



Charged pions (on right in figure [7]) can be produced in direct reaction at a rate larger than neutral pions only for 
T > 35 MeV. The photon channel dominates. 

IV. DISCUSSION AND CONCLUSIONS 

We found that the production of ir° is the dominant coupling of electromagnetic radiation to heavy (hadronic) 
particles with m 3> T, and as we have here demonstrated that noticeable particle yields can be expected already 
at modest temperatures T € [3,10] MeV. In present day environment of 0.1 -1 J plasma lasting a few fs, our 
results suggest that we can expect integrated over space-time evolution of the EP 3 fireball a tt° yield at the limit 
of detectability. For T — >• 15 MeV the tt° production rate remains dominant and indeed very large, reaching the 
production rate R' ~ 10 15 [MJ _1 fs -1 ]. Charge exchange reactions convert some of the neutral pions into charged 
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FIG. 7: On left: Muon and on right charged pion production rates in electromagnetic processes normalized by n production 
rate. Solid line (blue) for 77 , dashed line (green) for e + e~ induced process. 



pions which are more easy to detect. 

In this situation it is realistic to consider the possibility of forming a chemically equilibrated fireball with ir° , , ^ 
in chemical abundance equilibrium. The heavy particles are produced in early stages when temperature reached 
is highest. Their abundance in the fireball follows the fireball expansion and cooling till their freeze-out, that is 
decoupling of population equation production rates. The particle yields are than given by the freeze-out conditions, 
specifically the chemical freeze-out temperature Tf and volume Vf, rather than the integral over the rate of production. 
In this situation the heavy particle yields become diagnostic tools of the freeze-out conditions, with the mechanisms 
of their formation being less accessible. However, one can avoid this condition by appropriate staging of fireball 
properties. 

The present study has not covered, especially for low temperature range all the possible mechanisms, and we 
addressed some of these issues in the introduction. Here we note further that the production of heavy particles 
requires energies of the magnitude to/2 and thus is due to collisions involving the (relatively speaking) far tails of a 
thermal particle distribution. If these tails fall off as a power law, instead of the Boltzmann exponential decay p^ |. 
a much greater yield of heavy particles could ensue. There could further be present a collective amplification to the 
production process e.g. by residual matter flows, capable to enhance the low temperature yields, or by collective 
plasma oscillations and inhomogeneities. 

These are just some examples of many reasons to hope and expect a greater particle yield than we computed here in 
microscopic and controllable two particle reaction approach. This consideration, and our encouraging 'conventional' 
results suggest that the study of ir° formation in QED plasma is of considerable intrinsic interest. Our results provide 
a lower limit for rate of particle production and when folded with models of EP 3 fireball formation and evolution, 
final yield. 

It is of some interest to note that the study of pions in QED plasma allows exploration of pion properties in 
electromagnetic medium. Specifically, recall that 1.2% fraction of ir° — > e + e~7 decays, which implies that the 
associated processes such as e + + e" — > "f + ito are important. We cannot evaluate this process at present as it involves 
significant challenges in understanding of ttq off-mass shell 'anomalous' couping to two photons. 
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The experimental environment we considered here should allow a detailed study of the properties of pions (and 
also muons) in a thermal background. There is considerable fundamental interest in the study of pion properties and 
specifically pion mass splitting in QED plasma at temperature T <L Am and in presence of electromagnetic fields. We 
already have shown that due to quantum statistics effects, the effective in medium decay width of 7r° differs from the 
free space value, see figure [5] In addition, modification of mass and decay width due to ambient medium influence 
on the pion internal structure is to be expected. Further we hope that the study of pions in the EP 3 fireball will 
contribute to the better understanding of the relatively large difference in mass between 7r°and tt^ . The relatively 
large size of the PE 3 environment should make such changes, albeit small, measurable. 

The experimental study of tt° in QED plasma environment is not an easy task. Normally, one would think that 
the study of the ir° decay into two 67.5 MeV 7 (+ thermal Doppler shift motion) produces a characteristic signature. 
However, the ir° decay is in time and also in location overlapping with the plasma formation and disintegration. The 
debris of the plasma, reaches any detection system at practically the same time instance as does the 67.5 MeV 7. The 
large amount of available radiation will disable the detectors. On the other hand we realize that the hard thermal 
component of the plasma, which leads to the production of ir° in the early fireball stage, is most attenuated by plasma 
dynamical expansion. Thus it seems possible to plan for the detection of ir° e.g. in a heavily shielded detection 
system. 

The decay time of charged pions being 26 ns, and that of charged muons being 2.2 /zs it is possible to separate in 
time the plasma debris from the decay signal of these particles. Clearly, these heavy charged particles can be detected 
with much greater ease, also considering that the decay product of interest is charged. For this reason, we also have 
in depth considered all channels of production of charged pions and muons. Noting that practically all charged pions 
turn into muons, we have also compared the production rates of 7r° with all heavy particles, see dot-dashed (green) 
line in figure [71 This comparison suggests that for plasmas at a temperature reaching T > 10 MeV the production 
of final state muons will most probably be by far easier to detect. On the other hand for T < 5 MeV it would seem 
that the yield difference in favor of 7r° outweighs the detection system/efficiency loss considerations. Future work 
addressing non-conventional processes will show at how low T we can still expect observable heavy particle yields. 

An effort to detect ir° directly is justified since we can learn about the properties of the plasma (lifespan, volume 
and temperature in early stages) e.g. from a comparative study of the ir° and ^ production. We have found that 
at about T > 16 MeV, the pion charge exchange 7r°7r° — > ir + ir~ reaction for chemically equilibrated 7r° yield is faster 
than the natural ir° decay, and the chemical equilibration time constant, see the dot-dashed line in figure [5] Thus 
beyond this temperature the yield of charged pions can be expected to be in/near chemical equilibrium for a plasma 
which lives at, or above this temperature, for longer than 100 as. 

In such an environment the yield of ir° is expected to be near chemical equilibrium, since the decay rate is compen- 
sated by the production rate, and, within 100 as, the chemical equilibrium yield is attained. Moreover, the thermal 
speed of produced n can be obtained from the nonrelativistic relation ^m(v 2 ) = |T, thus v oc VT and, for T = 10 
MeV, v ~ 0.5c. This is nearly equal to the sound velocity of EP 3 , v s ~ cj\/2> = 0.58c. Thus the heavy n° particles 
can be seen as co- moving with the expanding/exploding EP 3 , which completes the argument to justify their transient 
chemical equilibrium yield in this condition. 

The global production yield of neutral and charged pions should thus allow the study of volume and temperature 
history of the QED plasma. More specifically, since with decreasing temperature, for T < 16 MeV, there is a rapid 
increase of the relaxation time for the charge exchange process, there is a rather rapid drop of the charged pion yield 
below chemical equilibrium — we note that charge exchange equilibration time at T = 10 MeV is a factor 10 5 longer. 
We note that the study of two pion correlations provides an independent measure of the source properties (HBT 
measurement). 

The relaxation time of electromagnetic production of muon pairs wins over 7r° relaxation time for T > 22 MeV, see 
dashed line, red, in figure O the direct electromagnetic processes of charged pion production (thin green, solid line for 
77 — > 7r + 7r~ and dashed, blue for e + e~ — * tt + tt~) remain sub-dominant. Thus for T > 22 MeV we expect, following 
the same chain of arguments for muons as above for charged pions, a near chemical equilibrium yield. If the study of 
all these tt°, tt^ , /x ± yields, their spectra and even pion correlations were possible, considerable insight into e _ ,e + ,7 
plasma (EP 3 ) plasma formation and dynamics at T < 25 MeV can be achieved. 
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